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O ! Abstract 

CN ■ We explore the relationship between the Glasma framework and the String Percola- 

\ tion Model by introducing a simple model for the rapidity distributions in nucleus- 

nucleus collisions. The model is solved for both symmetrical and asymmetrical col- 
•■O . lisions. The phenomenon of limiting fragmentation is briefly discussed. 
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O ' 1 Introduction 

The String Percolation Model (SPM) pQ for high energy hadronic collisions has 
a long history of successfully describing experimental observations [2f3Hi5] . As 
with any efficacious phenomenological model, one would wish, if not to fully 
derive it from first principles, to comprehensively understand its relationship 
rS [ with more fundamentally motivated approaches. 



The Glasma, a highly coherent configuration of classical longitudinal chromo- 
electric and chromo-magnetic fields, has been proposed [6f7f8f9p0lll2llll] as 
an adequate description of the matter formed immediately after a high energy 
hadronic collision. Further, it may provide |7.8.9 lT2fTT] the necessary physi- 
cal mechanism for the rapid thermalization observed at RHIC and thus for 
the eventual creation of the Quark Gluon Plasma. The Glasma hypothesis 
follows from the Colour Glass Condensate (CGC) approach to high energy 
hadronic processes [T3|14lll5lll6|17|18|19|20|21l 22.23.24.25j. which in turn pro- 
vides a theoretically sound and phenomenologicaily successful description for 
the infinite momentum wavefunction of a hadronic object. The CGC/Glasma 
framework tackles the same physical issues as those addressed by the SPM. 
The relationship between the CGC/Glasma and the SPM has been hinted at 
repeatedly, but it is only recently that efforts have been made to understand 
and explore such a connection [T2l26f27j . 
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In the CGC/Glasma framework, a nucleus-nucleus collision, viewed in the 
reference frame of the centre of mass of the collision, is to be thought of 
as that of two Lorentz contracted to the transverse plane sheets of coloured 
glass [28 29 30,31 32,33,34]. The CGC description of the infinite momentum 
hadronic wavefunction relies on the introduction of a longitudinal momentum 
separation scale. The small-x component of the wavefunction, a system of 
non-perturbatively large colour fields, can be described in terms of a classical 
Weizsacker- Williams field radiated by the large- x components of the wavefunc- 
tion. Such an effective description, in which the relevant dynamical degrees of 
freedom ( small- x gluons) are described by a classical field with the large- x con- 
stituents playing the role of a classical static (in light-cone time) source, holds 
for modes with longitudinal momentum close to the separation scale. The sep- 
aration scale is, however, arbitrary and can be lowered in successive steps to- 
wards increasingly smaller x. This renormalization group procedure is encoded 
in the B-JIMWLK evolution equations [Tfi|17lll8lll9f2Uf21ll22|l2~3 24 35 36 3 T1I3"8"] 
for the CGC. The small- x part of the wavefunction is characterized by a 
transverse momentum scale — the saturation momentum Q s — arising from 
finite density induced non-linear interactions between the wee partons (glu- 
ons). Modes with transverse momentum below Q s are saturated and cannot 
be further populated. Any additional gluon emission will thus result in an 
increase of the saturation momentum. 

Immediately after the collision, the hadron-hadron system is described by the 
highly coherent state of the two CGC sheets. The original transverse fields, 
associated with the large-x components of the hadronic wavefunctions, remain 
intact and, as the sheets pass through each other, add to themselves a distri- 
bution of colour electric and colour magnetic charge. These charge densities, 
which are of equal magnitude but opposite sign, act as sources for longitudinal 
fields living in the region between the parting CGC sheets. 

The physical situation described by the Glasma, the system of purely longi- 
tudinal fields in the region between the parting hadrons, is analogous to that 
underlying the SPM. As the effective number of strings, including percolation 
effects, is directly related to the rapidity particle density, we thus have the 
'equations' 

[Glasma] ~ [Effective Strings] ~ [Particle density] . 

It can thus be argued that once the Glasma is created through the mechanism 
outlined above and carefully detailed in [6lfTTj , its evolution can be phenomeno- 
logically described in the framework of the SPM. In this paper we follow such 
a line of thought and address the question of the rapidity evolution of a system 
of longitudinal fields/strings. 

In Sec. [2] we introduce our simple model and compute the density distribution 
in rapidity for both symmetrical and asymmetrical nucleus-nucleus collisions. 
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In Sec. |3] we address the issue of limiting fragmentation. Our conclusions are 
presented in Sec. HI 



2 Rapidity distributions 



2.1 The model 



The simplest model that can be imagined to describe the generation of lower 
centre of mass rapidity objects from higher rapidity ones with asymptotic 
saturation is the nineteenth century logistic equation for the dynamics of pop- 
ulations [391H0] 

1 :(P-V), (1) 



d(-A) 5 

where p = p(A, Y) is the particle density, Y is the beam rapidity, and 

A=|y|-Y. (2) 

The variable — A plays the role of evolution time. The parameter 5 controls the 
low density p evolution and must therefore depend on intrinsic parameters of 
the theory. The parameter A, responsible for saturation, will naturally depend 
on the nature of the target and on the energy e Y . 

The Y-dependent limiting value of p, determined by the saturation condition 



dp 



8(-A) 



0, is given by 



Py = \- (3) 



The condition qtzm2 = 0, defining the separation between the region A > 
A of low density and positive curvature and the region A < A of high 
density and negative curvature, gives 

Po = p(A ,Y) = ^. (4) 



Integrating JTJ between p and some p(A), and using (J3J) and (jlj), we obtain 

p(A,Y)= A _ff . (5) 

Eq. (jSj), illustrated in Fig.HJ is nothing but a generalization of the Fermi 
distribution, and is known to fit RHIC data for rapidity distributions in the 
case of central nucleus-nucleus colisions [4"l~ll4"2"] . There are three parameters in 
(jSJ): two from the original differential equation, py and 5, and an additional 
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one, Ao, from the integration of ([I]). The parameters py and A are expected 
to be energy (or Y) dependent, but S can be taken, at least approximately, as 
a constant. In the real Fermi distribution one has py = 1, but here we are not 
dealing with a dominantly fermionic system. 

The separation, at A , between a dilute region (A > A ) of objects with 
rapidity close to the beam rapidity and a dense region (A < A ) of low 
rapidity objects is a close analogue of that underlying the CGC framework. 
There, one distinguishes between a dilute set of objects (partons) carrying a 
large fraction of the hadronic wave function total longitudinal momentum - 
the source — and a dense region of partons with small longitudinal momentum. 
A further dense-dilute separation, that set by the saturation momentum Q s , 
is considered in the CGC. Such a separation is absent from our simple model 
since all dependence on tranverse dimensions has been neglected. 




-y a o 

A=|y|-Y 

Fig. 1. The distribution ([5]) showing saturation at large —A, the change of curva- 
ture at the point Ao which separates between the dense (A < Ao) and the dilute 
(A > Aq) regions. 



2.2 Symmetrical nucleus-nucleus collision 

So far we have concentrated on following a single nucleon pair. In a nucleus- 
nucleus collision there is a number N part of nucleons that take part in the 
collision. In the symmetrical situation of a AA collision we have 

N part = 2N A , (6) 

where N A is the average number of participants from each nucleus. Both in 
the CGC approach and in the SPM the particle density is proportional to Na 
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and we can write 

p(A,Y,N A )= Zl PY ■ (7) 

In order to justify (jTJ) we shall look in more detail at the SPM. The particle 
density is directly related to the effective number of strings 

(i 71 

p(A) = Nff ~ F( V )N S (A) , (8) 

where N S (A) is the number of strings contributing in a given A, F(rj) is the 
colour summation reduction factor 143 



and r\ is the percolation transverse density parameter 

V =Q 2 N S (A). (10) 

Here, r is the transverse radius of the string and R is the radius of the over- 
lapping area which can be obtained from simple geometrical arguments 



r l/3 

where R p is the proton radius. 



r~r p nx\ (n; 



The percolation transition occurs for [H] r] = i] c ~ (1.15 — 1.5). 

We distinguish two regimes: 

i) High density regime [A < A , r] = t] NaNa > r] c ] 

In the dense regime physics is dominated by gluons and sea quarks and the 
number of strings contributing is proportional to the number of collisions 
that is 



Further, rj is large such that F{rf) — >• r\ l l 2 and (see ([9 



N s ~ NT . (12) 



^)~T^' ( 13 ) 

so that we obtain 

p(A)~N A . (14) 
ii) Low density regime [A > Aq, r] = 7] pp < r] c ] 

In the dilute region physics is dominated by the valence strings and thus 
N s ~ Na, F(rf) is independent of Na and we obtain [4"6|47f48| 

p(A)~N A . (15) 
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We find that, at least approximately, the relation p ~ Na holds in both limits 
thus justifying <^7§. 

2.3 Asymmetrical collision 

We consider next the asymmetrical situation of AB collisions. In what follows 
we shall always assume that A < B and that A moves from left to right. Let 
us start by making a simple estimate, based on geometrical arguments, of the 
number N co u of nucleon-nucleon collisions resulting from a given number of 
participant nucleons (see Fig.[2]). If a tube of N 1 / 3 nucleons from nucleus A 

1/3 

interacts with iV B nucleons from nucleus B along the collision axis then, as 
the overlapping area of interaction is N^ 3 , we have 

N coll ~ (N^N^NT = N A N l J 3 . (16) 

Clearly, the well known result from Glauber calculus N co u ~ A^ 3 is recovered 
for a symmetrical situation. 




Fig. 2. Estimate of number of collisions for an asymmetrical AB collision (see text). 

The quantities N 1 / 3 and N l J 3 depend on impact parameter: for b ~ 0, N 1 / 3 ~ 
A 1 / 3 and N l J 3 ~ B 1 / 3 ; for very large b, N 1 / 3 ~ N l J 3 ~ 1. In general, A"j!/ 3 and 

1/3 

Afg can be theoretically calculated or experimentally determined, by using 
forward/backward calorimeters to detect the non-interacting nucleons: (A — 

2/3 1/3 

Na) in the forward direction and (B — N A ) in the backward direction. 
We distinguish again two regimes: 

i) High density regime [A < A ] 

In this case all the strings (valence and sea) contribute, their numbers 
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being proportional to the number of collisions 

N s ~ Ncou- NaN 1 / , (17) 

and we obtain, from (jSJ), 

p(A ~ — y) ~ N% 6 Nb /6 . (18) 

ii) Low density regime [A > Ao] 

The number of strings is essentially proportional to the number of valence 
strings, which is in turn proportional to the number of participating - 
forward or backward — nucleons, 

N S (A ~ 0) ~ Np ar t (Forward/Backward) , (19) 

We obtain, see (J5J), 

p(A ~ 0) ~ N A (Forward) , (20) 

and 

p(A ~ 0) ~ N% 3 Nl /3 (Backward) . (21) 
Bearing in mind (J7j), (TT8"j) . (I2"0l and (l2Tj) we arrive at 

p(A) = / PY (Forward) , (22) 



e^s— + 1 



p(A) = ; * " ^ (Backward) . (23) 

Na\ 



^V /6 e^ + l 



If we impose, in (1221) and ([23"]) . the condition A^^ = iv"^ (i.e., a symmetric 
situation) we recover (0) and the forward (1221 and backward ( 1231) distributions 
become mirror images. Let us see, in this simple situation, how the forward and 
backward distributions can be matched. The variable A starts, for y = ±Y, 
at zero and decreases towards the mid-rapidity region reaching a minumum 
value. We thus have the constraint 

- A B MIN - A F MIN = 2Y . (24) 

In the symmetric case Af IIN = A^ IIN and 

A MIN = -Y, (25) 

corresponding to y = and the distribution is centred around the centre of 
mass of the nucleon-nucleon system. 
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In the asymmetric case the denominators in ([22]) and ([23]) are different, and the 
condition (12^1) is not sufficient to match the distributions. Further requiring 

A^--ln^ = A^ + -ln^, (26) 

we obtain, after combining ([24]) and (!26l) . 

-Af //iV = r + ^ln^, (27) 

-A* JAr = Y-|ln^. (28) 
6 ATa 

The conclusion is simple: the maximum of the distribution p(A) is displaced 
towards the backward region (higher nucleonic density) and for the same A, 
Ps(A) > pa (A). These facts are qualitatively seen in the data 



We shall now be more specific regarding the quantities 5, py and Ao- We 
expect 5 to not have a strong dependence on Y, being essentially a constant. 
The parameter py is the normalized particle density at mid-rapidity and thus 
should be related to the gluon distribution at small x, 

Py ^ e XY , (29) 

with A ~ 0.25 — 0.3 [501I5T] . From energy-momentum conservation A was also 
estimated: A ~ 2/7 — 1/3 [12|52] . Regarding A we cannot assume it to be 
constant, as in the Feynman- Wilson gas model. As py increases with rapidity 
Y (1291) . energy conservation arguments give that A has to decrease linearly 
with Y. We write 

A = -aY , (30) 

with < a < 1. 

Eq. ([70 can be rewritten as 

1 , A , 2 dn e XY , . 

— p A, Y) = —— = A+aY . (31) 
N A N part dy e —i^ + l 

In fact, ( 13T1) was used in [H] as a suitable parametrization of RHIC Au — Au 
data from [531. 



3 Limiting fragmentation 



Recently, the phenomenon of limiting fragmentation [51] has received renewed 
attention [12,53, 55,56,57). The particle density p is, in general, a function 
of both A and the beam rapidity Y . The limiting fragmentation hypothesis 
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essentially states that for A larger than some F-dependent threshold, the 
density p becomes a function of A only. Since A is a decreasing function of 
Y (1301) . limiting fragmentation can be recast, in the framework of our model, 

as 

Pi^Y)y IZTf( A )- ( 32 ) 

Y -+oo , A>Ao 

From (13TT) we can write 

p(A, Y) — -> e -^ e (x- a /s)Y ? (33) 

Y -^oo , A>Ao 

so that limiting fragmentation will hold if, for 5 > 0, one has 



a 

-r 1 - (34) 



RHIC data [53] does not seem to favour (|34|) as recently observed in 
the tendency being for the curve, as the energy increases, to approach a step 
function. However, the BRAHMS collaboration [SS] points out that the net 
baryon rapidity loss seems to saturate between SPS and RHIC energies, which 
certainly is an indication for the limiting fragmentation regime. In our opinion, 
no strong statements should be made as the considered range in Y is small 
and QCD corrections most certainly introduce some weak violation of limiting 
fragmentation. 

If we assume (1511) to hold, that is that limiting fragmentation should hold in 
some limit, we arrive at the conclusion that the evolution of p(A, Y) with Y 
is governed by the same evolution equation as its evolution with A (pQ). 

In fact, from (131]) and (134j) . one easily arrives at 

1 e- A ' 5 

+ (35) 



which can be obtained from (1311) through the change 

XY < y-A/5. (36) 

Eq. (|36|) means that the evolution equation in Y is, see (JT]) and (j3J, 

dp(A,Y) 



\(p-—pA, (37) 
V Pa J 



BY V p A 

where pa = e~ A//<5 (obtained from py ( 1291) by using ( 1361) ) is the maximum 
value of p. If one rescales p to a function with maximum value 1, i.e. perform 
the transformation p — > pa p, then (1371) can be rewritten as 

^ = *(!-,). (-) 
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which is formally equivalent to the Balitsky-Kovchegov (BK) equation [351159] 
for the Y evolution of the scattering amplitude of a colour dipole once the 
dependence on transverse coordinates has been neglected. 



4 Conclusions 

We explored the connection between the Glasma and an effective description 
in terms of the SPM by introducing a simple model for the distribution in 
centre of mass rapidity of the particle density. 

The solution of this model for a symmetrical nucleus-nucleus collision (j7j), i.e. a 
collision of identical nuclei, reproduces (1311) the form used as a parametrization 
of RHIC Au - Au data [H]. 

In the case of an asymmetrical collision the maximum of the distribution 
is, in qualitative agreement with RHIC data , displaced towards backward 
(higher nucleonic density) region. The expression obtained for the forward (1221) 
and backward (123]) distributions, supplemented by the matching conditions 
( 1271) and (1281) . can be used, since it also reduces to (|3T|) for the limiting case 
of a symmetrical collision, to perform a comprehensive fit of RHIC data and 
to yield predictions for Pb — Pb collisions at the LHC [SD] • 

Our model yields, in terms of its parameters, a necessary condition (I3"4"|) for 
the occurrence of limiting fragmentation. This condition is not favoured by 
current RHIC data. The fit of RHIC data performed in [12] yields values for 
the parameters (A = 0.247, a = 0.269 and 5 = 0.67) that do not fulfil (l34j) . 
Further, attempting to perform a fit with (I3"4"|) as a constraint results in a poor 
description of the data. It should, however, be noted that, in the fragmenta- 
tion region, (|3~T|) with parameter values obtained without enforcing limiting 
fragmentation (as above) results in curves |42j that are indistinguishable from 
a universal behaviour. Our model predicts a clear violation of this apparent 
universal behaviour in the fragmentation region for Pb — Pb collisions at the 
LHC. 

To the best of our knowledge, (134)) first written in |42j and obtained here from 
a physically motivated explicit model, is the only possible indication that the 
apparent universal behaviour of the fragmentation region in RHIC data may 
not correspond to limiting fragmentation. 

If limiting fragmentation is enforced into our model, that is the fulfilment of 
( IMl) is imposed, than an equation (137|) for the evolution of the particle density 
with total rapidity Y can be written. This equation can, in turn, be rescaled to 
an equation (I3"5j) formally identical to the BK equation with no dependence on 
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transverse degrees of freedom. A note of caution as to the scope of this anal- 
ogy is due. The equation obtained here (|38|) applies to the string/gluon field 
configuration existing immediately after the collision and describes the evolu- 
tion with Y of the multiplicity of strings (i.e., particle density) of any given 
rapidity length A. CGC inspired studies of limiting fragmentation [56|57] rely 
on that gluon production in the fragmentation region of each of the colliding 
nuclei can be described, at leading order, in the framework of /^-factorization 
such that the inclusive cross section for gluon production can be written as 
the convolution of the unintegrated gluon distributions of the two nuclei. The 
kinematical region of fragmentation is such that the relevant Bjorken-x do- 
mains of the fragmenting nucleus (projectile) and the target are very different: 
large Bjorken-x for the projectile and small Bjorken-x for the target. The BK 
equation is then used to evolve the target to small-x (large Y), while the 
large-x gluon distribution of the projectile is obtained from a phenomenologi- 
cal extrapolation |57j . The BK dynamics guarantees that the typical transverse 
momentum in the projectile is much smaller than that in the target and that 
unitarity of the target evolution is preserved. Limiting fragmentation follows 
naturally from the target reaching its black disk limit and the gluon produc- 
tion cross section being determined by the universal large-a; component of the 
projectile. 

From the above discussion it should be clear that our formal analogue of BK 
(1551) and the BK equation as used in CGC based studies of limiting fragmen- 
tation apply to different stages of the collision and should not be compared 
directly. If, however, in future works one finds an evolution equation with 
Y for particle density distributions which is formally analogous to the BK 
equation with no transverse dependence, then this will trivially lead to limit- 
ing fragmentation. Whether such a strong statement can be made for formal 
analogues of the full BK equation or for analogues of equations beyond the 
mean field approximation of BK is not clear to us and clearly merits further 
consideration. 
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